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The Ginzburg-Landau free energy functional for hard-spheres is constructed using the Fundamen- 
tal Measure Theory approach to Density Functional Theory as a starting point. The functional is 
used to study the liquid-fcc solid planer interface and the properties of small solid clusters nucleating 
within a liquid. The surface tension for planer interfaces agrees well with simulation and it is found 
that the properties of the solid clusters are consistent with classical nucleation theory. 

PACS numbers: 64.10.-|-h,68.08.-p 

I. INTRODUCTION 

The nature of the liquid-solid interface is of both practical and theoretical significance. From a practical point of 
view, the interfacial surface tension is, along with the bulk free energies, one of the determining factors controlling 
nucleation of crystals from solution and plays as well an important role in the closely related phenomena of wetting. 
However, a complete description of the liquid-solid interface based on fundamental considerations remains a challenging 
theoretical problem since it requires knowledge of the spatial -dependence of the free energy in highly inhomogeneous 
systems. Much can be learned from more phenomenological approaches such as phase-field theory and model free 
energy functionals but a detailed relation of interfacial properties to molecular interaction models requires more 
fundamental approaches. On the other hand, very detailed approaches such as molecular dynamics and Monte Carlo 
computer simulation can provide direct microscopic information about interfaces, but theory is required in order to 
understand the results obtained. 

Classical density functional theory (DFT) is based on the fact that the Helmholtz free energy is a unique functional 
of the average density profile [3, It has been used to study liquid-solid interfaces - both coexistence and wetting 
- almost since its inception jJl^^Q. The most commonly studied system is that of hard-spheres since many density 
functional theories work best in this case. Although somewhat artificial, the hard-sphere interaction plays an important 
role in equilibrium statistical mechanics since more realistic pair interactions can be described based on perturbation 
theory about the hard-sphere interaction^ or by developments within DFT inspired by perturbation theory j^. 
However, the theories which most accurately described the bulk liquid and solid phases for hard-spheres suffer from 
a serious technical deficiency in that they allow for configurations in which the spheres overlap - a situation which is 
excluded on physical grounds and which should manifest itself as a divergence in the free energy for configurations 
in which the spheres touch f6l|. Thus, the liquid-solid interface can only be studied with these theories if either the 
allowed densities are artificially restricted, as in iXs ^] , or if ad hoc modifications are made to the DFT so as to create 
the required divergence at overlap0,0. In recent years, a new class of DFT models, generally known as Fundamental 
Measure Theory (FMT), have proven successful in a number of applications particularly involving inhomogeneous 
fluids (e.g., fluids near walls) [l0|,[lll, [13. The FMT approach has the advantage that, by the nature of the models 
used, the problem of overlapping hard-spheres is automatically solved - overlapping spheres lead to an infinite free 
energy penalty as one would expect. The goal of this paper is therefore to revisit the problem of the liquid-solid 
interface for hard-spheres both to further test the generality of the FMT approach and also as a preliminary step 
towards the study of interfaces for more realistic systems. 

Within the DFT-FMT framework, there is still considerable latitude in the level of description of the physical 
system. A bulk liquid is characterized by a constant average density p (T^) — p while a bulk solid is characterized by 
a spatially-varying density which must have the symmetry of the underlying crystal lattice. Inhomogeneous systems 
with either interfaces between different phases or confined geometries necessarily have more complex density profiles 
involving non-periodic spatial variations. For liquids, this does not not pose too great a challenge but for solids, the 
superposition of the spatial variations within a unit cell and the larger-scale variations coming from the interface can 
lead to numerically- intensive calculations see e.g ref. 0- Indeed, a recent application of the DFT-FMT approach 
to the problem of hard-sphere liquid-solid coexistence has been made and serves to illustrate the difficulty of this 
approach .13). In this work, the goal is to use a reduced description whereby a small set of order parameters rather 
than the fully detailed local density 14]. In this sense, the present work follows the spirit of phase field methods. 
However, the free energy functional is systematically derived from the microscopic description thus eliminating the 
need for phenomenological assumptions and thus allowing for quantitative predictions. 

In the next Section, the elements of DFT and its relation to the Ginzburg-Landau theory are reviewed. The FMT 
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density functionals are also described and the use of these functionals to calculate the elements of the GL free energy 
functional is presented. In Section III, the properties of the planer liquid-solid interface are determined both by means 
of parametrized profiles of the density and crystallinity and by numerical solution of the Euler-Lagrange equations. 
The calculated surface tension is shown to be in reasonable agreement with the results of molecular dynamics and 
Monte Carlo simulations. The structure and free energy of small solid clusters is also discussed. It is found that 
the properties of the critical nucleus - its size and the widths of the interfacial region - are well predicted by the 
results obtained for the planer interface using classical nucleation theory. The last Section summarizes the results and 
discusses possible refinements of the calculations. 



II. THEORY 



A. Density functional theory 

Density functional theory is based on the fact that the grand potential, the thermodynamic free energy appropriate 
for a system with constant chemical potential /i, constant temperature T and fixed applied one-body field can 
be written as 

n^F[p]~ J pp{~r)d~r + j <j) (T^) p ) dl^ (1) 

where the first term on the right is a unique functional of the local density. Different applied fields will give rise to 
different density profiles and the fundamental theorem of DFT says that there is a one to one correspondence between 
applied fields and density profilesj3|,0. It follows from this definition that the Hclmholtz free energy. A, is 

A^F[p]+ J (j>{~r)p {~r) d~r (2) 

so that F [p] is the intrinsic contribution to the Hclmholtz free energy due to the density profile. For situations in 
which the applied field is not important - for example, when it only represents the walls of a container, the second 
term is unimportant in the thermodynamic limit and F [p] is often referred to simply as the Hclmholtz free energy. 

At fixed field, temperature and chemical potential, the density must minimize the grand potential giving the Euler- 
Lagrange equation 

F[p]-p. + ^{r)^Q. (3) 
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Given a model for the intrinsic free energy functional, this gives a closed description of the system which 

generally takes the form of an integral equation. Most applications of DFT make use of parametrized density profiles 
so as to reduce the effort needed to determine the density. In general, if the density is given by p (~r*) ~ p (7^; F), where 
/o(T*;r) is some fixed function of the spatial coordinates and the permeates F ~ {ra}"^li , then the Euler-Lagrange 
equations become 

where the notation indicates that F[p\ = F (F) is an ordinary function of the parameters. The derivatives are 
understood to be evaluated holding all parameters constant except Fa- Note that if one of the parameters corresponds 
to the average density, say Fq = p, and if the field is zero (or confined to the boundaries so that it can be neglected 
in the thermodynamic limit) this gives 

^F[p] = p (5) 
dp 

^/[P] = 0, a>0. 

The first equation is just the usual relation between the Hclmholtz free energy and the chemical potential while the 
second says that the Helmholtz free energy must be stationary with respect to all of the other parameters. 

For a bulk solid corresponding to a Bravais lattice with a single atom per unit cell, the density must have the 
symmetry of the lattice and so takes the form 
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where | i?„| are the lattice vectors. In this case, the average density is 

P = j P{~r)d~r (7) 

= Piatt / f{'r )d~r 
Jws 

where the second integral is restricted to the Wigner-Seitz cell and where "Piatt is the lattice density, defined as the 
number of lattice points per unit volume. The integral therefore defines the occupancy, = p/Pzatt- occupancy 
of one means that every lattice site is occupied, a value less than one means that there are some vacancies, a value 
greater than one means that there are some interstitials. Note that eq.® can equivalently be written in terms of 
Fourier components as 

p(r) = ^exp(ii?„.r)/p„) (8) 

n 

where f (^k^ is the Fourier transform of / (V) and where |_ftr„| is the set of reciprocal lattice vectors. 

A typical parametrization of the density, widely used in practical calculations, is to take / {~r^) to be a Gaussian so 
that 

/(y) = Xo(^)'^'exp(-ar2). (9) 

There are then three parameters: the average density, p, the width of the Gaussian, a, and the lattice density JJiatt 
and the occupancy is just xq — J)/ piatt- In this case, one has / {j^ri^ = pcxp (— _ft'^/4a) which, together with the 

fact that Ko = shows that lim^^o P {~^\ F) = p so that the parametrization of the density encompasses both the 
Gaussian approximation for the solid and the uniform liquid. For this reason, it is common to take the value of the 
first non-trivial Fourier component to be a measure of the " crystallinity" , denoted m, giving 

771 = exp {-Kl/4,a) (10) 

so that 771 = corresponds to a uniform fiuid and 77i = 1 to an infinitely localized solid. Real solids have values of 77i 
which are close to, but always less than, one. 

To study interfacial properties, it is necessary to allow for spatial variation in both the average density and in the 
crystallinity. Here, we follow Ohensorge, et al|a| and allow p and a to vary giving 

p (y ) = (p ) ipi^u) ('^\ ' f-" (^) - (11) 



or 



p (y ) = p (y ) (^Kn ■ exp {-Kll^a )) (12) 

n 

depending on which form of the bulk density is used as a basis for the generalization. These expressions are obviously 
not equivalent although it will turn out below that within the Ginzburg-Landau framework, the differences are 
unimportant. An alternative parametrization used by Haymet and Oxtoby[3,0| is 

p (^ ) = 5p (r ) + j2 f-« (^) - Rn) ') (13) 



or 



p (^) = Sp ) + ft,,, Y exp {^Kn ■ exp {-Kl/4a (^ )) (14) 

n 

where it is assumed that p^ — 'Pi^n in the solid. One difficulty with this form is that since Sp (V) must be allowed to 
be negative (the liquid is less dense than the solid) there is the unphysical possibility that p (V) < for some points 
"r*. The form given in ea. (|ll|l is positive definite thus avoiding this problem. 
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B. Fundamental Measure Theory 

The original form of Fundamental Measure Theory as given by Rosenfeld[l3,[lll is in some sense a development of 
scaled particle theory. This was further extended by Tarazona, Rosenfeld and others using the important requirement 
that the known, exact free energy functional be recovered in the one-dimensional limit of the theory [H, [13, d. 
The resulting class of theories has proven successful at describing inhomogeneous hard-sphere fluids including the 
hard-sphere solid. The theory involves a number of local variables, tIq which are linear functionals of the local 

density of the form 

na{~r i) = j d~r 2 p{'r' 2) Wa [1^12) ■ (15) 

The set of weights Wa (~?^12) include a simple step function 8 (f- — r) which serves to define a local packing fraction 

V (^1) = y rf^2 P (T^2) e (I - na) (16) 

as evidenced by the fact that in the uniform limit, p{~r2) ^ P, one has 77 (~r^i) = ^'pcr^ which is the usual definition 
of the hard-sphere packing fraction. All of the remaining weighting functions are tensorial dyadics of the form 
rr...r5 (r — ^) and the resulting variables are written generically as 

Tij...i (Vi) = J d~r2 ri2,iri2j...ri2,iS (^ri2 - p(^2) • (17) 

It will be useful to give simpler names for the first two of these quantities, namely 

sCri) = j d~r 2 5 {r 12 - p {1^2) (18) 



Vi{~ri) = j d~r 2^12, iS [r 12-'^ p{~r 2) 



where the names stand for "scalar" and "vector" respectively. 
The Helmholtz free energy functional is written as 

F [p] - F,a [p] + [p] (19) 

where the ideal part of the free energy is 

md [p] ^ JdV ip (V) Inp (V) - p (V)) (20) 
with f3 = l/iksT), and the excess contribution is written in the FMT as the integral of a function of the local variables 



[p] ^ ({n„ (T)}) (21) 



which is usually expressed as 



with 



= </'l + 02 + 03 (22) 

= -^s(r)ln(l-7y(r)) (23) 
TTcr- 

^ 1 s"^ {~r) ~ {~r) 

The form of ^3 depends on the particular version of FMT. Here, three common theories will be considered. The 
first theory was proposed by Rosenfeld et aL12j and is perhaps the simplest form of FMT capable of giving a good 
description of the hard-sphere solid 



7r(l-?7(^))' V s^ r) 



■ (24) 



5 



The second is the theory of TarazonaT^ which makes use of a tensor variable 

/3<A^ = ^ ^^_,^^2 (^) • (y ) . ^ (r ) - s ) (y ) - Tr (7^3 )) + 5 (r ) Tr {j^ )) ) . (25) 

Both of these theories have the property that they reduce to the Percus-Yevik approximation for the hquid. The third 
theory builds in the more accurate Carnahan-Starling equation of state via a heuristic modification of the Tarazona 
theory^^iES- It is commonly known as the "White Bear" functional and is given by 



_ 2 7y(7) + (l-r7(7))-ln(l-ry(r)) ^^ 
*3 - 3 ^f<^) 



(26) 



All three of these theories have deficiencies. The RLST theory incorporates the Percus-Yevik approximation for the 
liquid which is not very accurate at the density of liquid-solid coexistence. To the extent that the RLST theory 
gives a good description of liquid-solid coexistence (see below), it is because it gets the liquid and solid "equally 
wrong" . The Tarazona theory also reduces to the Percus-Yevik approximation for the homogeneous fluid but it gives 
a better description of the properties of the homogeneous solid leading to a poor description of coexistence|l7|. For 
this reason, this theory has not been used in the present investigation. The White Bear (WB) functional gives an 
improved description of the dense fluid by incorporating the Carnahan-Starling equation of state in an ad hoc fashion. 
As a result, the implied pair distribution function for the fluid, obtained via the Ornstein-Zernike equation, will not 
vanish in the core region as it should. Nevertheless, the free energy does diverge for overlapping hard-spheres as in the 
other forms of FMT. The conclusion is that the RLST is probably the best theory in terms of the formal properties 
of the free energy while the WB may be expected to be the better in terms of quantitative results. 



C. Ginzburg-Landau theory 

If it can be assumed that the order parameters vary slowly over atomic length scales, a simplified free energy 
functional can be rigorously derived from the exact free energy functional by means of a gradient expansion in the 
order parameters. When carried out to second order, the result takes the form of a phcnomcnological Ginzburg-Landau 
free energy functional[3,|l,|i,[l3,|23,[ll, 



13^0 L [F] = dR 



(r (r)) - p^Tp (r {r)) + (f {r 



dVa {r) {r) 



dR^ 



dR, 



(27) 



where the mean field term PF (F) is the free energy of a bulk system evaluated with order parameters F and 



The coefficient of the gradient term is 



1 

2V 



p(F)= / p(ri;F)d1 



d~r'id~r'2 7-i2i'"i2jC2 (Vi, ^^2; P) 



dp{-t^-V) 9p(V2;P) 



dVa 



dVb 



(28) 



where the direct correlation function is also evaluated for a bulk system with constant order parameters and is 
determined from the free energy via the standard relation 



C2 {'ri,~r2\Y) 



S^I3F,, [p] 



5p{r i]T)5p{r 2;F) 



(29) 



As discussed in ref.(14|. this expression is derived under the assumption that the lattice structure is held fixed. This 
precludes the use of any characteristics of the lattice, such as the lattice constants or the primitive lattice vectors, as 
order parameters. Nevertheless, it is possible to explore the liquid-solid interface since, as discussed above, this can 
be done for a fixed lattice structure. 

When the underlying lattice has cubic symmetry, as will be the case here, and in a coordinate system aligned with 
the principle axes of the lattice structure, the symmetry under 90 degree rotations around the axes implies that the 
coefficient of the gradient term can be written as 



K 



ah 



gabSij + hab {Six^jy + SiySjj; + Si^^Sjz + SizSjj; + SiySjz + SizSjy) 



(30) 
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where 



9ab (r) = ^ J dr idr2 ^1202 ( r 1, r 2; T) — (31) 

1 ^,dp{r,;T)dp{r2;r) 



hab(r) = ^ y"d^lC?^2 ?'l2x?'l2yC2 (^1,7*2; r) 



■ 

To apply this formahsm to the FMT, note that the direct correlation function takes the form 

C2[r 1, r 2; r) = - 2^ / dr — Wa [r - r i)wi3(r - r 2) (32) 

At first, it appears that the evaluation of the matrix Kf^ will involve triple volume integrals making it extremely 
expensive to calculate. However, it is in fact possible to arrange the calculate so that it takes no more effort to 
calculate Kf^ than is required to calculate the free energy. To do this, note that the calculation requires evaluation 
of the integral 

.->.-> (-^ -> \ ^ .dp{~ri]T)dp{~t2]T) 

d r id r 2 ri2iri2jWa [ r - r x)wp,{r - r 2) 7:^ 7^ (33) 

Ci a Ci b 

which can be simplified by introducing some additional functionals of the density. Specifically, let 

(7*) = y dT^i riiWa (7^ - 7^1) p (T^i; r) (34) 
i^a C^) = j d'r'irurijWai'r -l^i) p{~fi;T) 



so that the full expression for Kf^ becomes 



i^ab L V dV ^'■^ " + ^r: arr^ (or^ 



2V^J dn^dnp (t>) dn^g (t>) 3< (?) dn}, (t>) 

This is of the same structural form as the expression for the excess free energy, namely a spatial integral involving 
functions of the local density functionals, and can therefore be evaluated as easily. Furthermore, note that for the 
tensorial densities, one has 

Tt.iC^) = f J) T,,,.,„, (y) (36) 
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™" (^) = (I) r,,..,™„(y) 



The only other quantities needed to evaluate the GL free energy functional are 77' (T*) and rj^^ C^)- Thus, for the 
RLST theory, the tensorial quantities must be evaluated up to third order while the White Bear theory requires the 
fourth order tensor as well. Explicit expressions for all of these quantities are given in Appendix^ An interesting 
analytic result, proven in Appendix IbI is that in the RLST theory hab (F) — 0. Using the WB theory, it is possible 
that hab (r) but the present calculations for the case of an FCC lattice are that the calculated values are so small 
that neglect of hab (r) makes no difference to the results reported below. Thus, for all practical purposes, the GL free 
energy functional is not affected by the anisotropy of the underlying lattice for hard spheres. 



D. Numerical methods 



Many quantities, such as the density, can be evaluated either in real space or in Fourier space. Generally, the 
former is more efficient in the (large a) crystalline state and the latter in the (small a) liquid state. For the present 
study, both methods have been used in the calculation of the density and of all of the reqiured density functionals (see 
Appendix^ for details). The implementation was checked by comparing both methods for values of aa^ ^ 20, where 

1 /3 

a — {4:/ Piatt) is the FCC lattice paramter. In all subsequent calculations, the Fourier-space method was used for 
aa^ < 20 and the real space method otherwise. 
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The spatial integrals were evaluated in a limited volume. Specifically, in a homogeneous system, one has that 

1 
V 



dT/30({n„(r)}) = / dr/30({n„(r)})/ / dr 

■11 J cell 



where the integrals on the right are restricted to the conventional unit cell. For the FCC lattice, elementary symmetry 
considerations show that the integral on the right can be restricted to the volume x, y, z > and x > y with a symmetry 
factor of 16. Further restrictions are possiblej^J but were not used. The integrals were then evaluated using an evenly 
spaced grid of 20 points in each direction ( e.g x — nSx/21 for < n < 20) . Increasing the number of points had no 
significant effect on the calculations. 

These calculations are still time consuming, especially at intermediate values of a, so both PF (F) and K"^' (F) 
were evaluated over a grid of points in parameter space and bi-cubic spline interpolation[22 used in the subsequent 
calculations. Note that since all of the elementary density functionals are linear in the average density, it is possible 
to perform the calculations for fixed a but many values of the average density at once, i.e. in parallel, thus saving 
time. One problem is that these quantities, especially Kfj' {T), are divergent at high average densities since, for 
sufficiently high average density, the local packing fraction rj (^)can exceed one which is just the signal of overlapping 
hard-spheres. For a ho mogeneous crystal, the value of the local density for which this divergence occurs is given by 
Tiatt/v where ri(r) is defined in Appendix 1X1 and given explicitly in ea. (jA6|) . This maximum density is always 
greater than p/^itt and approaches in the limit a ~^ oo. Since most of the variation as a function of p occurs near 
the maximum value, the procedure used was to discretize the density as 



Pk 



Pmax 2P™^^'"^P fcma"'"' ) ) ' ^ fcmax/2 



with pniax chosen to be slightly below Piatt/vi^)- the calculations reported below, 200 points were used for the 
density and the crystallinity was sampled on an evenly spaced grid of 80 points in the range < m < 0.95. 
Some analytic checks on the numerical calculations are possible. Note that the density can be written as 

P (^) P^'i^P exp (iKj ■ ^) (37) 

j 

= mi^V^^^) exp {iKj ■ 

3 

= p + pm ^ exp (^iKj ■ +pto(''2/''i) ^ exp (^iKj ■ T"^ + .... 

j&Si ]eSs 

where Sn is the set of n-th shell reciprocal lattice vectors. It immediately follows that, in the liquid limit, one has 

m^O dp 

y 9p(y;F) f.-y 

lim = P > exp Uil,- 

m^o dm < ^ ^ V ■' 

giving 



OC 



lim Opp(F) - \i, I r^c',{r;p)dr (39) 

lim 5p,„ (F) = 

T 2 _2.^ f°° i I / sinfcir 

liniffmm(r) = -^T^p Ni I r C2{r;p)— dr. 

m^o 3 Jq kir 

where C2 (ri2;p) is the direct correlation function of the liquid, ki is the magnitude of the radius of the first shell 
of the reciprocal lattice and A^i is the number of eleme nts i n the first shell. For an FCC lattice in real space, the 
reciprocal lattice is BCC giving iVi = 8 and Ki — (^) With the RLST theory, the dcf of the liquid is just the 

Percus-Yevik dcf. 



(r^2--p) = + af ^ + af (^) j 9 (r - a) , (40) 
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with coefScients 

\2 



PV ^ 

^2 = 2'^O'0 



2 



whereas the White Bear functional gives the same functional form but with coefficients 

_ [ 2 - rj + Utj^ - 6ii^ ^ 21n(l-7?) ^ 

_ f 3 + 5v{v-2){l-rj) , 21n(l-r;) ^ 

Equations H39|l - H43() give a useful check on the full numerical calculation. 

Except for the cases explicitly discussed above, all numerical integrals, minimizations and the solution of ordinary 
differential equations were performed using routines from the Gnu Scientific Library • One-dimensionsal minimiza- 
tions were perfomed using either Brent's method or bisection while numerical integrals were performed using adaptive 
integration (the GSL "QAGS" routine'23'|) with relative and absolute accuracies set to 10"'*. Multidimensional mini- 
mizations were performed using the Simplex algorithm of Nelder and Mead, see e.g. ref. as implemented in GSL, 
which was terminated whent the simplex size was smaller than 10^^. 



III. RESULTS 



A. Bulk coexistence 



As a baseline for the interfacial calculations, it is necessary to know the predictions of the various theories for the 
liquid-solid transition. Using the Gaussian parametrization for the density, eqs.®-©, the expressions given above for 
the free energy are evaluated and, in accord with eq.(|SJl, minimized with respect to a and the lattice density pjatt while 
the average density must be adjusted so as to satisfy the relation between the free energy and the chemical potential 
given in eq.©. One minimum is always found at a = corresponding to the uniform liquid having density (p). 
Another occurs at ag > corresponding to a bulk solid having some average density (/i). The true equilibrium 
is whichever of these solutions that minimizes the grand potential: bulk coexistence occurs when they give identical 
values 

^pF {-pi ifi) , 0) - p^ipi (m) = ^pfiF (p, i^i) , a.) - (3^lp, ifi) (44) 

or, using eq.® and the definition of the pressure, /3P = 

f3P{Pi{fi),0)^(3P{pAf^),as) (45) 

which is the usual condition of coexistence. The results for the RLST and WB theories are shown in Table^as well as 
the values from previous calculations and the values from simulation. These numbers are on the whole constant with 
those in the literature, particularly when it is noted that the calculations reported in the literature were performed 
with the occupancy fixed to be one. The one exception to this general agreement is the RLST theory where the 
present results differ noticeably from those reported by Rosenfeld et al. 12]. However, as an independent check, the 
evaluations of Warshavsky and Songj^^ are also shown and are seen to be consistent with the present calculations 
thus supporting the accuracy of the present evaluations. 
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TABLE I: The solid density 'p^giO* , liquid density 'PugCr'^, reduced pressure pPa^ and chemical potential ^ at bulk coexistence 
as determined from (a) the present work, (b)Rosenfeld et al [l^ . (c) Warshavsky and Song [l^ . (d) Roth et al.|ljj| and (e) 
from the simulations of Hoover et al|24ll. 



Theory 










RLST(a) 


1.017 


0.935 


12.18 


16.88 


RLST(b) 


1.031 


0.938 


12.3 


17.05 


RLST(c) 


1.020 


0.937 


12.26 


16.99 


WB(a) 


1.022 


0.932 


11.20 


15.67 


WB(c) 


1.023 


0.934 


11.29 




WB(d) 


1.023 


0.934 






MD(e) 


1.040 


0.940 


11.70 





B. The planer interface 



The simplest application of the theory is to the planer liquid-solid interface. The order parameters are assumed to 
vary in only one direction, which is writen as u = ^ ■ n where n is the normal to the interface, from a bulk solid in 
the limit u — oo to a bulk liquid at u ^ oo. For these calculations, the GL free energy functional takes the form 



(3n 



GL 



/oo 
du 



1 



P^^P + 2 (9ab (r (u)) + 2hab (r (u)) {n^ny 



nxTiz + riyn^ j) T'^ (m) T'^ (u) 



(46) 

where A is the area in the x-y plane. The calculations were performed for a fixed lattice structure at the lattice 
density appropriate for bulk phase coexistence as determined above. 

It is important to realize that the liquid-solid interface is only stable for a unique value of chemical potential. This 
is because the interface is expected to be a localized structure that decays towards the bulk phases exponentially 
as one moves away from the interface|25). However, since the bulk free energy function, f]F{T), and the matrices 
gab (r) and hab (r) are being approximated via interpolation of tabulated values, it is unlikely that the coexistence 
properties will be exactly the same as found above. Therefore, a necessary step is to determine coexistence at this 
fixed lattice structure using the interpolated f3F (T). The procedure used was as follows. Given an initial guess of 
'Psoi, the value of asoi was determined by finding the minimum of (3F (j)g^i,Q). From this, the chemical potential is 
determined from the thermodynamic relation, eq.© where it is important that the derivative is taken at constant 
a and p/att- Then, the liquid density is determined from the requirement that /9F(p;j^,0) give the same chemical 
potential. Finally, the bulk pressures, -^/Si^ (F) — are compared for the liquid and solid and the difference used 
to adjust the value of Pj,^;. The process is iterated until a value of p^^; is found that results in equal pressures between 
the bulk phases. The results of the bulk calculations are given in Table ITTl where the "exact" calculations and those 
based on the interpolaton scheme are seen to be in good agreement. 

One interesting quantity displayed in this table is the difference between the lattice density and the average density. 
The values given for the RLST theory are in reasonable agreement with the asymptotic values obtained analytically by 
Groh 26J. However, the WB theory gives a small but negative difference in densities which would mean that instead 
of vacancies, the theory predicts interstitials. This surprising and unphysical result could be an artifact of the WB 
theory (which, as discussed above, is an ad hoc extension of Tarazona's theory) or it could be a numerical artifact 
indicating that the true vacancy concentration is too small to be reliably calculated with the numerical methods used. 
Decreasing the tolerances of the numerical evaluations in these calculations did not result in a change of the sign of the 
vacancy concentration leaving open the possibility that the result is a real artifact of the theory. Only analytic work 
along the lines of ref.|2^ will resolve this question. Nevertheless, it is important to note that the chemical potential is 
very sensitive to the differences between p^^; and 'piatt i or other words the value of the occupancy. This is because 
the chemical potential is given by 



d 1 



(3F{V) 



d 1 



(3F{T) 



+ 



d 1 



Psol d^latt V 



PF{T) 



(47) 



Working with a fixed value of pj^^j means that the second term on the right is neglected. This is not a problem if the 
free energy is stationary with respect to the lattice density as it should be at thermodynamic equilibrium. However, 
for typical values of a in the solid, the maximum average density, defined as that at which the free energy and all 
other quantities diverge, is only on the order of 0.0001% above the lattice density so that the free energy is very 
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TABLE II: The solid density p^^iU'^, liquid density PnqCr"^, difference between lattice and average densities, reduced pressure 
pPa^ and chemical potential p. at bulk coexistence as determined from the "exact" calculations and from the numerical tables 
with interpolation. 



Theory P,^;^" Ph,'^" Piat^ " P.oi'^' M PPo'' 

RLST-exact 1.017 0.935 4 x lO"'^ 16.88 12.18 

RLST-interpolated 1.017 0.937 3 x lO"'^ 16.98 12.27 

WB-exact 1.022 0.932 -7 x 10"** 15.67 11.20 

WB-interpolated 1.021 0.934 -2 x 10"** 15.75 11.28 




FIG. 1: The crystallinity (a), and the average density (b), profiles as functions of position perpendicular to the interface as 
determined using the RLST DFT and parameterized profiles. Shown are results using hyperbolic tangent profiles with Bp = Bm 
(solid lines), allowing Bp ^ Bm (dashed lines), and with a Gaussian term in the density profile (dotted lines). 



sensitive to changes in the average density and the lattice density. Thus, neglecting the difference between these and 
setting p^q; = p/^jt, can lead to errors on the order of 10% in the chemical potential when evaluated at fixed lattice 
density even though the free energy itself is insensitive to this small difference. 

The structure of the liquid-solid interface is determined by minimizing the free energy functional, ea. H46l) . with 
respect to the spatially-dependent order parameters leading to the Euler-Lagrange equations 

^-b (r (u)) — Tb {u) + — — — — -BF (r {u)) ~p^ip\=Q (48) 



du^ ^ ' V 2 dTa J du du dTa\V 

with 

Kab (r) = gab (r (m)) + 2hab (F {u)) {nxTiy + n^n^ + riyri^) 

and the boundary conditions 

hm r, - (49) 

u — > — oo 

hm Fa = ri 



where = {nisoi, Psoi): etc. Direct solution of the Euler-Lagrange equations is difficult due to the presence of 
unstable solutions. So, in addition to this method (discussed below), parameterized forms of the order parameters 
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FIG. 2; Same as Fig. [Hbut using the White Bear DFT. 
were explored. Perhaps the most natural choice is to model the order parameters as simple hyperbolic tangents, 

r,(u) = + (r: - ri) i (tanh {a^ {u - + 1) (50) 

so that one minimizes the free energy with respect to the widths and positions of the interface. The surface tension 
(actually, the surface free-energy) is calculated as 

r = 1 [pn^^l--^ [T] - pn'^l') (51) 

where ^^q^ is the free energy of either bulk phase. The results are shown in Table imi First, the simple hyperbolic 
tangent profiles with the interface fixed at i3„i = Bp — Q gives surface tension of 0.73 for RLST and 0.75 for WB. 
Allowing the interfaces to move gives a relative displacement of nearly one hard-sphere diameter and lowers the surface 
tensions to 0.67 and 0.66 respectively while giving a considerably narrower density profile. As a further refinement, 
inspired by the numerical results discussed below, a Gaussian term was added giving 

Va=T\ + {Tl - ri) i (tanh {A^ (u - BJ) + 1) + exp (u - E^f) . (52) 

The addition of the Gaussian makes not appreciable difference for the profile of the crystallinity, but can affect the 
profile of the average density. The profiles of the order parameters are shown for the RLST and WB theories in 
Figs. mi and Q respectively. It is clear from the figures that there is little change in the profile of the crystallinity 
but the density is sensitive to functional forms used. These results were obtained using as initial values Aa = O.Scr, 
Ba — Ca ~ Ea ^ and Da — 1 . Variation of the initial values by a factor of two has relatively little effect on the final 
result. However,, starting with much sharper widths, Aa > 2.5a, leads in some cases to a different minimum with a 
slightly lower surface tension and in other cases no minimum is found. It could be argued, however, that these cases 
are unphysical as the sharpness of the variation of the density violate the assumptions underlying the derivation of the 
Ginzburg-Landau form of the free energy (namely, that the order parameters vary slowly over atomic distances) . The 
planer density profile, p (u) ^ -j J P (^) S {u — 'r^ ■ n) (Pr , corresponding to the RLST result using the combination 
of hyperbolic tangent and Gaussian is shown in Fig.||2J)) for a profile in the [100] direction. It can be seen that in 
agreement with the results reported by Warshavsky and Song[l3, the interface involves perhaps 8 lattice planes. The 
planer density corresponding to the offset hyperbolic tangents is indistinguishable from the one shown. In fact, from 
Figs. Q and ||2Jl it is clear that while all three profile functions give indistinguishable results for the crystallinity, the 
offset hyperbolic tangents and the hyperbolic tangent plus Gaussian give very similar profiles for the density and differ 
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FIG. 3: The atomic density averaged over planes perpendicular to the interface as a function of position, calculated using the 
RLST theory and the offset hyperbolic tangent parameterization. The position is shown in units of the interplaner spaceing 
for [100] planes, d = 0.5a where a is the lattice parameter. 

from the simple hyperbolic tangent profile. Furthermore, while in all cases one finds that going from the bulk liquid 
to the bulk solid, ordering (i.e. an increase in crystallinity) precedes densification (an increase in average density) due 
to the signficantly greater width of the crystallinity profile, the more complex parameterizations accentuate this by 
shifting the densification curve towards the bulk solid region. 

The direct integrations of the Euler-Lagrange equations is numerically challenging as there are both stable,physical 
solutions and unstable, unphysical solutions and the accumulation of numerical errors means that eventually all 
numerical integrations of the equations become unstable. To see this, note that far from the interfaces, one expects 
that the Euler-Lagrange equations can be linearized about the bulk values of the order parameters giving 



where Fa are the bulk values of the order parameters and ST a = Ta — Ta . The possible solutions are determined by 
the eigenvalues and eigenvectors of the matrix occurring in this equation. It turns out that in both the solid and liquid 
regions, there is one positive eigenvalue, allowing for decaying solutions, and one negative eigenvalue, corresponding 
to undamped oscillations. It is the mixing of the unphysical oscillatory solution with the decaying solution that causes 
numerical difficulties as illustrated in Fig. Q for the RLST theory. The figure shows the result of integrating both 
from the bulk liquid towards the interface and from the bulk solid towards the interface with initial conditions chosen 
to allow for matching the solution at some point near the interface (thus constructing a shooting- method solution). 
Starting from the liquid, it is in fact possible to integrate almost completely through the interface until, in the solid 
region, pollution from the oscillatory solution causes the density to take on unphysical values. (Of course, unphysical 
values lie very close to the bulk value of the average density thus requiring relatively little inaccuracy to achieve this.) 
Note that the curve so obtained shows a slight depletion of the density near the interface as well as some structure 
on the solid side of the interface. Both of these effects are very small, given the overall small change in density 
between the liquid and solid, and are probably of little physical consequence. Integrating from the solid side shows 
similar effects as well as similar oscillations in the bulk liquid region. The surface tension of the shooting solution is 
"ya^ /ksT ~ 0.67 which is not much different from that of the curve beginning from the liquid or solid sides which give 
0.64 and 0.65 respectively. These figures are consistent with the surface tension obtained using the parameterized 
profiles and support the contention that the parameterizations are reasonably accurate. 

The surface tension of the planer interface has been determined from both Molecular Dynamics (MD) 27] and Monte 



Carlo (MC) simulations "^i]. The MD gives values of ja'^/ksT = 0.62,0.65 and 0.58 for the [100], [110] and [111] 




(53) 
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FIG. 4: The crystallinity, left, and the average density, right, obtained by numerical integration of the Euler-Lagrange equations 
based on the RLST DFT. In the case of the crystallinity, the hyperbolic tangent profile discussed above is also shown. For the 
average density, the curves obtained by integrating from both the bulk solid and the bulk liquid regions is shown. 



TABLE III: The order parameter profile parameters obtained by minimizing the free energy. The profiles studied are the 
hyperbolic tangents with Bm = Bp (H), the "offset" hyperbolic tangents where Bm 7^ Brho (OH), and the hyperbolic tangents 
with a Gaussian term (HG).Also included are the results from MD simulations of ref and the MC simulations of ref. 281 . 
In all cases, the last column gives the surface tension. 



Theory 


Profile 


-Am 


A, 


Bp 




Dp 




"fo^ /ksT 


RLST 


H 


0.61 


0.83 






* 


* 


0.730 


RLST 


OH 


0.67 


1.64 


-0.70 




* 


* 


0.669 


RLST 


HG 


0.68 


0.99 


* 


-0.039 


1.27 


0.04 


0.667 


WB 


H 


0.74 


0.84 




* 


* 


* 


0.754 


WB 


OH 


0.85 


2.54 


-0.78 




* 


* 


0.659 


WB 


HG 


0.88 


1.70 




-0.06 


1.97 


-0.21 


0.656 



MD 0.617 
MC 0.623 



directions respectively while the MC gives ^a^ /ksT = 0.64,0.62 and 0.61. These results show a weak dependence on 
the orientation of the interface and the observed variation is not consistent between the two methods which perhaps 
indicates that the asymmetry is very small. The agreement between these values and those calculated from the 
GL model are therefore quite good. As expected, the WB model, which incorporates a better equation of state for 
the liquid than does the RLST model, seems to be slightly more accurate. The remaining differences between the 
calculations and the simulations can at least in part be attributed to the imposition of an invariant lattice structure. 



C. Solid Clusters 



The GL theory can also be used to study the structure of solid clusters embedded in the liquid, 
symmetric system, the grand potential becomes 



For a spherically 



(F {R)) - P^rp (R) + \g., (F (i?)) 



R'dR 



(54) 
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giving the Euler-Lagrange equations 



_2 d „2 dTb (R) , fdgat {T (R)) 1 dgtc {T (R)) \ dn (R) dV, (R) 8 



(55) 

The boundary conditions are 

lini = (56) 

R^O dR ^ ' 

Um r„ (i?) = Y\ {R) 

it— »(X) 

Note that the first condition ensures that the first derivative along any Hne passing through the origin is continuous, 
e.g. 

Un.^:^.lin.^. (57) 

Classical nucleation theory (CNT) is based on the observation that, to a first approximation, the excess free energy 
of a solid cluster of radius R embedded in a large volume of fluid will be 

47r 

j^^Qcluster ^ ^cluster jpj _ ^liquid ^ _^ ^3^^b«;fe _^ 47ri?^7 (58) 

3 

where Aw*"*"^ is the difference in the bulk solid and liquid free energies per unit volume. The stability of a cluster of a 

given radius clearly depends on the value of ^Afl'^'"'**'^''. If the liquid is the favored state, so that U.q^^'^ > 17^'^*'', then 
J_^^c;«ster ^ Q fQj. g^jj ^ meaning that any cluster will collapse to i? = 0. If the solid is favored, /317g'|;"^ < jSn'^''^"^, 
clusters smaller than Rc ~ 2j (Aw''""^) ^ will collapse while those with R > Rc will be unstable towards unlimited 
expansion. At this point, ^s-Ail'^'"**®'' = — 87r7 < so the critical cluster is a maximum in the free energy with 

3 (Aujb-"'f 

This has been derived under the assumption that the surface tension is independent of R. In fact, as defined above, 
the surface tension only applies to conditions of coexistence so that for hard spheres it is a unique property. Another 
way of expressing this is that this classical nucleation model can only apply very close to coexistence. In that case, 
one can also write 

A^bulk ^ (^p^ _ _ (^^p^ _ ^^^^ (go) 

) + {Ps- 

n + ipi - pt^ 




'''-l^)Ps) 
= (p,-p;)(/,— -^) 




\ I ^coex\ ,.coex i I Z? t' ^coex\ i ^coex\ ,.coex 

) + [Ps- Ps )P - MPs - -rji'l [Pi ) + (Pl- Pi )P - m 



giving the well-known expression 



3 {Ps - Pl) {p. 



As in the previous Section, a parametrized form for the order parameters is considered. It is again assumed that a 
hyperbolic tangent is a reasonable guess for the shape of the interface so, taking account of the boundary conditions 
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FIG. 5: The excess free energy verses the lattice packing fraction, rjiat 
the WB theory (open symbols) for Rp = 30 and (HAfi = 0.25 



T^piattcr'^ /& for the RLST theory (filled symbols) and 



suggests using 



(Fa (0) 
(Fa (0) 



1 + baR 1 - ta-nh {AaiR - Ra)) 
l + {baRf 1 - tanh (-Aai?a) 

l+baR ( 1 + exp {-2AaRa) 



(62) 



1 + {h^Rf V 1 + exp (2Aa (R - Ra)) 



with 



ba = 



2A„ 



1 



(63) 



This function has vanishing gradient at i? = and for large R decays as exp {—AaR) / R which is the expected 
asymptotic forni[2^. Note that the radius of the cluster is determined by Rp and i?m which are of course not 
necessarily equal . For a fixed value of the chemical potential, /z, the values of Fjj are computed from the known 
properties of the bulk liquid. Then, the profile is used in ea. (|54|l and the free energy minimized with respect to 
Ta (0) , A,n,Ap and Rm for a fixed value of Rp, which is taken to define the cluster size. In order to avoid unphysical 

and rp(0) = pmax - 



regions, such as to > 1, auxiliary variables u and v are defined via r,„(0) = rrim 
where rrimax and pmax are the maximum values occurring in the tables. 

The critical cluster then corresponds to the value of Rp for which the free energy is stationary. One difference from 
the planer interface is that the properties of the solid cluster cannot be assumed to be the same as those of a bulk 
solid with the specified chemical potential. Since the cluster is of finite size, the order parameters at the origin, Ta (0), 
will in general differ from those of the bulk solid and, most importantly, the lattice parameter cannot be assumed 
to be that of the bulk solid. Thus, in addition to minimizing with respect to the parameters of the profile, it is also 
necessary to minimize with respect to the lattice parameter. These points are illustrated by Fig[Sl showing the excess 
free energy as a function of the lattice density for particular values of the cluster radius and supersaturation. As 
expected, the excess free energy shows a minimum for a particular lattice density. Figure El shows the variation of 
the resulting lattice densities as a function of the radius of the cluster. For small clusters, the density is significantly 
lower than that of a bulk solid at the same chemical potential and increases rapidly as a function of cluster size. For 
larger clusters, the rate of change decreases as the bulk limit is approached although the variation with cluster size 
is still noticeable even for the largest clusters (i?i — 170cr). The WB theory gives slightly higher densities than the 
RLST theory, as would be expected from the bulk coexistence data (see Tabled). In all cases, the average density at 
the core, p{0), is very close to the lattice density. 
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FIG. 6; The lattice packing fraction, riiatt = T^piatt'y'^ /& as a function of the cluster size for the RLST theory (filled symbols) 
and the WB theory (open symbols)for /3Ap = 0.25 

Figures I3a) and I3b) show the excess free energy as a function of the cluster size for different values of the 
supersaturation determined using the RLST and WB theories respectively. Also shown are the predicted excess free 
energy based on CNT, eq. (|58|) . using the surface tension at coexistence obtained from the planer calculations of 
the previous Section. Clearly, CNT gives a reasonable approximation to the structure and energy of the critical 
cluster. In fact, this is a very sensitive test of the agreement between CNT and the detailed calculations. The largest 
discrepancy observed in the figures occurs for the lowest supersaturation, /3A/x = 0.125 and at largest cluster sizes. 
This is surprising since large clusters approach the bulk limit while surface tension becomes less important and one 
would expect that the CNT calculation would become increasingly accurate. However, since the radius of the critical 
cluster diverges as /3A/.t 0, the CNT results are very sensitive to the value of the chemical potential for small /3A/i. 
This is illustrated in Fig. [7Jb) where nearly perfect agreement between GL-DFT at /3A/x = 0.125 is found with the 
CNT result for /3A/.t = 0.120 thus suggesting that the differences seen are at least in part due to numerical inaccuracies 
in the determination of the chemical potentials at coexistence. Other factors contributing to the disagreement are 
that the CNT assumes that the cluster has the properties of the bulk solid and that the surface tension is that for 
bulk coexisting (/3A/i = 0) phases. 

Figure |S| shows that the lattice densities for small clusters are well below that of the bulk lattice density. A 
significant difference between the two DFT's is apparent in that the RLST theory gives a discontinuity in the lattice 
density as a function of cluster size while the WB theory does not. The discontinuity arises because the free energy 
as a function of lattice density (for fixed cluster size and supersaturation) calculated using the RLST theory has two 
minima. For small clusters, the low-density minimum dominates while for larger clusters, the high-density minimum 
has lowest free energy. Further calculations have verified that for the lowest supersaturations shown in Fig. |S1 the 
high-density minimum does dominate for sufficiently large clusters. Furthermore, the calculations also confirm that 
the high-density minimum is the one for which the density and crystallinity are closest to that of a bulk solid. 

Finally, in all cases, the structure of the liquid-solid interface is consistent with the results for the planer interface 
with similar values for the widths and relative displacements of the crystallinity and average density curves. 

IV. CONCLUSIONS 

The primary results of this paper are the formulation of the Ginzburg-Landau free energy functional for hard- 
spheres based on Fundamental Measure Theory and the use of this functional to study the properties of the liquid-solid 
interface and of small FCC solid clusters in solution. It was shown that the required elements of the GL free energy 
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FIG. 7: Excess free energy as a function of the cluster size as determined using (a) the RLST theory and (b) the WB theory. 
Curves are shown for /3A/^ — 0.125, upper curves, to f3Afj, = 1.00, lowest curves. The predicted excess free energies from CNT 
are also shown (broken lines). 



functional can be calculated reasonably efficiently using FMT provided that the family of functionals of the density 
are extended. The resulting functional was used to study liquid-solid coexistence with a planer interface as well as 
the structure of small solid clusters. 

For the planer interface between coexisting liquid and solid phases, the resulting surface tension is in reasonable 
agreement with simulation. No dependence of the surface tension on the lattice plane was found which is also consistent 
with simulation. Furthermore, the results using parameterized profiles and obtained via direct integration of the Euler- 
Lagrange equations were found to be consistent, notwithstanding the numerical difficulties of the latter procedure. 
This is an important point as, e.g., the early result obtained by Curtin using the WD A DFtTI, which preceded 
determination of the surface tension for hard-spheres by computer simulation and seemed to be in good agreement with 
the later simulations, was subsequently shown by Ohnesorge et al.0 to be spurious and due, apparantly, to the use of 
over-constrained profiles. That same parameterization was also used by Marr and Gast 8] and, with some modification, 
by Kyrlidis and Brown,23|. Recently, Warshavsky and Song, hereafter WS, performed similar calculations using FMT 
without the Ginzburg-Landau approximation. They observed greater differences from simulation than in the present 
work and non-negligable spatial asymmetry. The agreement between the present results and simulation may be 
fortuitous or the difference between the present results and those of WS may be due to other details in the later 
calculation such as the assumption of unit occupancy. In any case, the agreement between both FMT calculations 
and simulation is an improvement over the results based on older DFT's which typically give a value of the planer 
surface tension that is about half of that measured in simulation jM 1^ l29l| . 

The GL functional was also used to study the properties of small solid clusters in superdense solution. For the 
range of supersaturtions considered here, it was not found to be possible to stabilize clusters of radius less than about 
15 hard-sphere radii; the free energy difference from the liquid is found to be very small and no local minimum in 
the free energy could be found. The results for clusters that could be stabilized are consistent with the predictions 
of classical nucleation theory. One interesting observation was that the RLST theory predicts a discontinuity in the 
structure of small clusters: very small clusters have unexpectedly low lattice densities while larger clusters approach 
the properties of the bulk solid. The crossover point between the two structures increases as the supersaturation 
increases and, conversely, appears to diverge near coexistence. However, since the RLST is based on a Percus-Yevik 
description of the fluid, which is not accurate at such high densities, and since the WB theory does not show this 
effect, it seems likely to be an artifact. 

For both the planer interface and the clusters, it was found that moving from the bulk liquid towards the solid, one 
first observes ordering of the fiuid and then densification. This is interesting as it is the opposite of the predictions 
of recent studies of crystallization directly from a low density gasJJOj. For the latter case, it seems to always be more 
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FIG. 8; The lattice packing fraction as a function of the cluster size as determined using (a) the RLST theory and (b) the WB 
theory. Curves are shown for /3A/i — 0.125, upper curves, to /3A/i — 1.00, lowest curves. 



favorable to densify first, thus forming dense liquid droplets, and then to order. Of course, the main reason that 
this scenario is not observed here is that for hard-spheres it is only possible to study nucleation of the solid from an 
already-dense fluid - there is no equivalent of the low-density gas-solid transition. This is also undoubtedly one of the 
reasons that CNT seems to work so well for hard-spheres (i.e., because there is no critical point). 

Although the Ginzburg-Landau model is reasonably successful in the applications described here, there are problems 
which cannot be discounted. Most important is that the free energy is unstable with respect to very sharp interfaces. 
Such rapid variations in the order parameters are outside the scope of the GL model, which is based on a gradient 
expansion, and so have been avoided here by always starting the minimizations from relatively slowly varying profiles. 
Within reasonable bounds for the starting parameters, the resulting profiles are then relatively insensitive to the 
starting point. 

While of some intrinsic interest in testing the GL model and the density functional theory, these results are only 
intended as baselines for more interesting studies of realistic interaction models such as the Lennard-Jones potential. 
Since the only input to the GL model is a reasonable model for the bulk free energy and a reasonable model for the 
bulk direct correlation function, rather than a full blown density functional theory, it is expected that the extension 
of this work to other potentials will be relatively straightforward. 
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APPENDIX A: DETAILS OF THE CALCULATIONS 
1. The density functionals 



All of the linear density functionals are of the form 

[p]) = j d'r'i w; (7* - p (7*i) . 
As long as the density is written as a sum of basis functions at each lattice site, 



then the functionals can also be written as 



with 



■ (^;[p]) = E^(^-^»;w) 



n{'r';[p])= / dT^i w (T* - T*i) p (^i) . 



(Al) 



(A2) 



(A3) 



(A4) 



Before proceeding, note that in the case of tcnsorial quantities, Ty...fc Cr; [p]), the only vector available is r and the 
only tensor, aside from fr, is the unit tensor. Thus, it follows that 



Vi (~r*) = V (r) fi 
T,,{r) = A{r)S,j+Bir)rirj 
Tiji (^) = C(r) {ri5ji + rj5ii + f;(5y ) + D{r)rirjri 

Tijlm {'r) = E (r) {5ij6lm + kl^jm + km^jl) 

+F (r) {nrjSim + riXiSjm + nrmSji + VjriSim + rjrmSu + rirmkj) 
+G (r) riXjrirm 



(A5) 



The basic functionals for FMT are then found to be 

-a (r - I) 



Vi^) = l(er/(^/S(r+|))-er/(^/S(r-|)))-^5(7 



(A6) 



v{r) 



1 / acr^ 



2r V IT 



exp I —a 



crar 



5(7") 



A{r) 


= V (r) 

crar 


B{r) 


= s{r) V (r) 

crar 


C{r) 


= —B{r) 
ara 


D{r) 


= v{r)-5C{r) 


E{r) 


a or 


F{r) 


= —v(r)-5E(r) 
ara 


G(r) 


= B{r)- 7F (r) 



(A7) 



Finally, one needs the additional quantities 

7?* (7*) = H{r)ri 
if^ (T) = I (r) fifj + J (r) 5ij 



with 



H (r) = rr] (r) — -^v (r) 
I(r) = rH(r) - ^B(r) 



In Fourier space one has that 



3 

Ha {^^^ = j exp (^—iK ■ 7^2^ Wa (^2) • 

The tensorial quantities can be expressed as in real space but with the vector r replaced by A; so 
Vi (^K^ = V (r) % 
Tij (if) = A (k) Sij + B (k) kikj 
~i]i(K^ = C{k) (k.Sji +%Sa +kAj^ + D{k)kikjki 

'^ijlm(^K^ = E {k) {Sij6im + SiiSjm + StmSjl) 
+F{k) 

+G (A;) kikjkikm 
The scalar and vector functional are 

s{k) = wa'^jo (^y) 
Vik) = ^TTcr^ (jo ( ^^o- ) +j2 ilka 



T,, 



2 ka . I k(T ^ 



v,{k) = -ma -n l^-j) 



and the coefficients for the tensorial functionals are 



B{k) = -7:a'h 



2 J \ 2 J 

ka 
~2 



= ia nj-j, 



ka 
~2 



-(') = -r^'(-(T)+-(T 

0(t) = A34(^ 
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Finally, 

= H{k)ki (A14) 
rf^ (k^ = I{k)kikj + J{k)Sij 



with 



56 V V2 J V2 



2. The derivatives of the free energy 

The free energy is written in terms of the integral of 

{{na (^)}) = ^1 {{na (^)}) + ^2 (K (r)}) + </.3 ({n„ (r)}) 

and for the gradient term in the GL functional one needs the second derivatives of this function. Consider the first 
two terms which are the same for all of the DFTs, 

^i+^^ = -i^^l-(l-^) + if^- (A16) 

The second derivatives of these are 

1 5 1 5^ — f ^ 

9r?2 TTCT^ (1 _ Tva (1 -r)f 



as"' TTCT (1 — T]) 

9^ . ^ 1 -1 r 
(01 + 02) = -; zOij 



and the cross-derivatives are 



(9^ Ills 

^(01 + 02) = + — 3 (A18) 

Or]OS TTCr^ 1 — ?? TTCT (1 — Tj)"" 

5^ . ^ I -V^ 

Q Q (01 + 02) = — -; 3 

^2 



^^^^^(01+02) =0. 

The third term takes different forms depending on the theory. For the RLST theory 

1 „3 / ,,2\3 

iRLST _ ^ 5 I . V 



and the second derivatives are 



<^ I RLST J- I ^ V 





247r 


1 


s3 


47r (1 




1 


s 


47r(i 


- vf 


1 


s3 



1-T2 (A19) 



a7?2^^ 47r(l-„)n 



1 - - (A20) 
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and the cross derivatives are 



drjds 

drjdvi 
52 



iRLST 
P3 



iRLST 
'3 



lRLST 



47r (1 - rjY 

1 sv^ 

1 Vi 



1 - 



1 + 



1- 



1+3- 



dsdvi 47r(i_7y)2 
For the WB theory, first consider the simpler Tarazona theory which has 

3 'v ■'t -1? - sv"^ -Tr (1^3^ + sTr ("t ^ j 



so 



167r 



dvidvj 



dT,,dTir. 



and the cross-derivatives are 



92 



7' 



drjds 
drjdvi 
drjdTij 

dsdvj^ 



^1 



92 



dsdTi-, 



dTijdvi 



Now, in the WB theory one has 



F{rj) 
F"{n) 



6 





3 "^ij ~t~ 



2s5i. 



3 Ti'. 



s6i. 



IGtt 



8^ {l-riY 



3 -36jiTmi - STjiSii + 2sSij5r, 



167r 



-u^ + Tr 



3 T^aVa - SUj 3 (Tia - S(5io) t^a 



47r {I-jjY 47r (1 - 77)' 

3 



87r(l-r?)^ 

3 SiiVj+SjiVj 
16n (1 - ,,)2 



2 ?? + (l-r;)'ln(l-r;) 
3 r?2 

2 _ 277 - 2 (1 - 77) In (1 - V) 

3 ^3 

2 677 - r;2 (6 - 4r;) In (1 - rj) 

3 rf 



(A21) 



(A22) 



(A23) 



(A24) 



(A25) 
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Then, 



= FMtJ-^^"^ (A26) 



dTadTb 



+ 5a7)5br,F" {ri) cf)^ 

<I>V^ = iv) + ^ (1 - V) (Sar, + St,) F' (,7)) ff^^<l>^ + Sa,St,F" {rj) (A27) 



APPENDIX B: VANISHING OF hab FOR THE RLST THEORY 
Combining eas. (|^ and gives 

^''^ (r) = ^ y r id r . '^i^^-i^. ( r - r u;, ( r - . ,) —^^ (Bl) 

The goal here is to prove that this quantity vanishes in the RLST theory. The idea behind the proof is to make a 
change of variables in the integral whereby = {xi,yi,Zi) — > T*^ = {—Xi,yi,Zi) for z = 1,2. This clearly gives an 
overall sign change as well as affecting the arguments of the various functions occurring under the integral. Using the 
fundamental fact that for a uniform solid (i.e. spatially constant T) with a cubic lattice structure, the density has 
reflection symmetry so that p{r i;T) = p(r '^;T) , it follows that the low order density functionals occurring in the 
RLST theory have simple parity and it is therefore possible to show that, aside from the overall change of sign, the 
integrand is invariant. This proves that hab (L) — 0. 
To begin, 

= ^l^jdr,dr,dr x,,y,, ^ ( , _ , J ( , _ , ^) __ 

- ^l^jdr,dr,dr x,,y,, ^ ( r - r,)w,{r - ^ .) 

a,f3 

where we have also made the change of variable — {x, y, z) — > T*' = (—2^, y, z). The claim is that A = 0, where 

A = ^A„/3 

A„, = anlW ' 

7\ 7\ Wa[r - ri)wi3(r - r 2) 

anaOnp 

Now, the scalar weight functions depend only on the magnitude of the separation, e.g. Wrj {5~r") — Wri {\5~r'\) = 
Wjj (|(5T^|), so that they are invariant under the change of variables. The vector density diJ"?*'!) — 
{SV/ I) e (|(5T"| - f ) = {SV/ \SV\) e {\51^\ - § ) so that Wy^ (|<5^'|) has odd parity and the other compo- 
nents have even parity. From this and the reflection symmetry of p{~r';T), it follows that the density functionals 
77 (V) , s (~r*) and 1? (T^) have the corresponding parities (that is , all are even except Vx (^') = —Vx C^))- From the 

explicit expressions for given in the previous appendix, it is seen that all of these have even parity except for 

those in which one and only one of the derivatives is with respect to Vx ("^)- It immediately follows that 

Ass — ^rji] — Ag^ — A^s — ^v^v^ — ^^y^y — ^v^v. 
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Finally, , — #73 , » a — and ^ a — are all of odd parity (because they are proportional to v^) as are 

WsWv^jWjjWv^jWv^Wvy and Wy^Wy^ so that A^^^ = A^„^ = Ay^y^ = A.y^y__ = thus proving that A = and hence 
that hab (r) = in the RLST theory. The same is not true of the WB theory as terms such as TjcTd, and therefore 



do not have a definite parity under reflection. 
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